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Abstract 

We consider the cosmological application of Lee- Wick theory where a field has a higher deriva- 
tive kinetic operators. The higher derivative term can be eliminated by introducing a set of 
auxiliary fields. We investigate the cosmological evolutions of these fields as a candidate of dark 
energy. This model has the same structure as so called "quintom' model except the form of poten- 
tials and the sign of the slope of the potentials. This model can give the stable late time phantom 
dominated scaling solution (wde < —1) or tracking attractors (wde = 0) depending on the choice 
of the slopes of the potential. In order to be a viable dark energy candidate, the present energy 
density contrast of dark energy (J^de) should be close to an observed value (0.73) at the same 
time. However, a simple toy model of the theory can not satisfy both J^l ~ -1 and Q^J, = 0.73. 
If we include the self interaction term of Lee Wick field, then we are able to produce the observed 
values of those quantities. However, in this case we are not be able to have the stable solutions 
and we need to suffer from the fine tuning of its mass. 



1 Introduction 



The instability of the Higgs potential due to the quadratic radiative corrections has motivated 
many extensions beyond the Standard Model (SM). Recently, there is a new extension of the SM 
so called Lee- Wick Standard Model (LWSM) [1], based on the work of Lee and Wick [2, 3] where 
they show a finite theory of quantum electrodynamics (QED) by including the regulator term as a 
higher derivative version of QED. If we extend the idea to all the SM field with a higher derivative 
term, then we have a LW partner for the every SM field. Its application to gauge bosons [4], 
neutrino masses [5], Maxwell-Einstein theory [6], fiavor changing [7], and electroweak precision 
data [8, 9, 10] have been considered. 

There is an almost identical work with the above proposal to demonstrate a single scalar 
field model of dark energy with a higher derivative shown the equation of state (eos) crossing —1 
[11, 12]. Similar idea adapted to inflation model is also considered [13]. By including the only 
self interaction term of the phantom field (M^'i/;^/2), this model can provide ilj^g = 0.73 and 
~ —1.0. However, this kind of potential sufl'ers from the fine tuning problem. Including 
the higher derivatives of a scalar field is well match to one of dark energy models, so called 
"Quintom" model [14] when we introduce a set of auxiliary fields. There have been many works 
done for quintom models with various potentials [15, 16, 17, 18, 19, 20, 21]. The perturbation 
of quintom models and their effects on cosmological observations are also studied [22, 23]. Even 
though there are potential problems of this model like the violation of causality and the instability 
of ghost field, we will not consider these aspects in this work. 

The energy density of the Lee Wick (LW) field dominates at late time because its energy 
density increases (a;^ < —1) while those of the scalar field and the barotropic fluid decrease 
(lo^,uj^ > —1) as the universe evolves. However, the eos of the effective dark energy (i.e. the 
combination of a;^ and uj^) can cross —1 depending on the slopes of potentials of two fields. This 
dark energy might be able to preserve the tracking behavior of the scalar field requiring less fine 
tuning at early universe and reach to phantom attractors at late time. Only self interaction of 
two fields are considered in many literatures except [17, 18]. However, we naturally have the 
interaction between the scalar field and the LW field in LW theory. 

In the next section, we will investigate the existence and the stability of tracker solutions of 
LW dark energy for the specific choice of potential by using the analysis of autonomous system for 
physical quantities. We check the cosmological evolutions of the dark energy and the barotropic 
fluid in section 3. We reach our conclusions in section 4. 

2 Lee- Wick Dark Energy 

To illustrate the effect of Lee- Wick theory in cosmology, we consider a theory of one self interacting 
scalar field, a, with a higher derivative term. The Lagrangian density is 
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= Iv^aV^a - - , (2.1) 

where Dcr = g^'^V /^Vucr. We can rewrite the above Lagrangian by introducing an auxihary scalar 
field ijj as 

Ca-^p = ^V^aV^a - ^Da + ^M^' - V{a) , (2.2) 

where Oa = M'^ip can be obtained from the above Lagrangian (2.2) to recover the original La- 
grangian (2.1). Now if we define a = (j) — tp, then the Lagrangian (2.2) becomes 

= ^V^(/)V^0 - ^V^i/jV"^j + ^mV' - V{(P - i^) , (2.3) 

where wc use the integration by parts [11, 12]. In this form there arc two kinds of scalar field, a 
normal scalar field (f) and an Lee- Wick (LW) field which has the opposite sign of the kinetic term 
compared to the normal field. Stability of this theory was considered in the literature [1]. We will 
present a phase-space analysis of these scalar fields in a homogeneous and isotropic flat universe 
containing a barotropic fluid with its eos {coj = 7—1) where the pressure and the energy density 
of the fluid is related by p-,, = (7 — l)p-y. 

As a toy model, we ignore the mass term in tp and assume the interacting terms of two scalar 
fields as an exponential potential. 

V{(P - V) = Vint(<^, V') = V^_^ exp(-/«A^<^ - KX^i;) , (2.4) 

where k = s/SttG. As shown in the above Lagrangian (2.3) we can not help having the interaction 

term between two scalar fields in LW theory. The homogeneous scalar fields (f> and ^p in a spatially 
flat FEW cosmological model can be described by the fluids with effective energy density and the 
pressure 



PDE = -4>^ - + Vint{(f>,'(p) = P4, + P4, -Vinti(l),ip) , (2.5) 
PBE = - ^'0^ - Mnt(^,'0) =P</> + ^int(0,V') = '^DEPDE , (2.6) 

where we indicate the effective eos of the dark energy as wde and define the energy densities and 
pressures of a scalar field and a LW field as 

Pi, = 7r + ^mt(0, 1p), Pi, = TT - Mnt(9i>, ^) = '^V/'V' • (2-8) 
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We show the second equahty of dark energy density and pressure related to those of a scalar field 
and a LW field because we have the interaction between two fields and only the total energy- 
momentum will be conserved. Now the Pridemann equation becomes 

= Y (^^°^ '^^'^ ^ y i^'^'f' ^P^~ ^int(^!', ip) + Pij = y Per- , (2.9) 

where pcr is the critical energy density and we again explicitly represent pde in the second equality 
to show the differences in the evolution equations with and without the interaction term between 
two scalar fields. Then the effective eos of the dark energy a;DE is given by 

;>DE </>2-^2-2Vi,t(</>,V') 

a;DE = = : • (■^■-LUI 

PDE </.2-V'2 + 2yi,t(0,V') 

This effective eos can be rewritten by using equations (2.7) and (2.8) as 

u^^ = Ul+''±A±^\^ (2.11) 
3 y iiDE / 

where Ode = We can find the evolution equations of two scalar fields and the fluid in fiat 

per 

FRW models by using the Bianchi identity which are 



^ + 3i/0 + ^^^4fl^ = 0, (2.12) 

o(p 

^ + SHi.-^y^4P^ = 0> (2-13) 
dip 

p^ + 3Hjp^ = 0. (2.14) 

We can find that the energy density of an individual scalar field is not conserved while the total 
energy density is conserved 

p^ + SHil + u;^)p^ = ^^"f^'^h , (2.15) 

p^ + 3H{l+u;^)p^ = (2.16) 

PDE + 3iJ(l + a;DE)pDE = 0. (2.17) 

To study the critical point structure and the stability of the system, it is convenient to introduce 
the following dimensionless variables [24] 



X-ip 



7!' 

Ktp' 

7f' 



z = 



K^/Vint{(^>, V') 



V3H ■ 



(2.18) 
(2.19) 
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Table 1: The values and the existence condition for fixed points. The various tracking attractor 
solutions are shown for the specific choices of A<^ and A^. We define A^^j = — 



Case Xij, y z Existence Ode 7de 

K - 1 A11A^,A^ 1 2 

F 1 A11A^,A^ undefined 

SP ^ ^ ^e//<6 1 % 
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where a prime means a derivative with respect to the logarithm of the scale factor In a. With 
these variables (2.18 - 2.19) we can express the Pridemann equation (2.9) and its derivative as 



12 2 I 2 I 2 

= x^-x^ + y +z 



H' ( 



7 2\ q 7 , /7 i\ 2 (1 i\ 2 , 7 2 

2')=^ -2 + (2-^)^^-(2-V^^+2^ 



(2.20) 
(2.21) 



Also we can find the evolution equation of those variables (2.18 - 2.19) 




y' = -yyj^ + ^i>^^x^ + x^x^)j , (2.23) 

4 = -(^ + 3ja;^-^Av,/. (2.24) 

Fixed points at finite values of x<^,x^ and y in the phase space are correspond to solutions where 
the scalar field (s) has a constant equation of state. Critical points correspond to fixed points 
where x'^ = y' = x'^ = and give an expanding universe with a scale factor given by a cx where 

"=6x2-6x2+37^2- (2-25) 

In Table 1, we classify critical points and show their properties and values. We show fixed 
points as one two dimensional hyperbola K which is kinetic dominated solutions, a fixed point 
F corresponding to a barotropic fluid dominated solution, a fixed point SP which is dark energy 
dominated solution, and a fixed point SPF meaning a fluid dark energy dominated solution. K 
and F points exist for any value of A^ and A^. However, there are constraints for the existence of 
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Table 2: Eigenvalues and stability 



Case 


mi 


1712 


7713 


Stability 


K 

F 




3 

\- 


3(2 - 7) 

3 

\- 


3 — ^~^{X^x^ + XijjX^) 

3 
2 


unstable 
unstable 


SP 


-3 + ^ 




^e//-37 


stable 
ifA2^^<37 



SPF 



■1(2-7) 



-1(2-7) 



(2-7)A 



-1(2-7) 



87(37-A^„) 



stable 



if37<A^,,<|^ 



points and SPF. The point SP exists only -when Ag^j = — < 6. Thus, if the slope of a 
potential of a scalar field A<^ is too big compared to that of a LW field A<^, then a point SP will not 

be existed. In this point the effective equation of state of dark energy is given by o^de = — 1- 
Also the point SPF exists only for A^jj > 87 cases. In this case, a scalar field evolves much faster 
than a LW field and it tracks the background fluid. Thus, the effective equation of state of dark 
energy follows that of a barotropic fluid wde = ^-y- 

We study the stability of critical points by substituting linear perturbations of variables (2.18) 
- (2.19) into the evolution equations (2.21) - (4.6) and find three eigenvalues for the evolution 
equations of the first order in the perturbations. Stability of system requires the real part of all 
eigenvalues to be negative. 

K : These kinetic dominated solutions always exist for any value of A^ and A^. It gives the 
equation of state of dark energy as 1 like inflation. This yields H'/H = —3. The linearized system 
of (2.21) - (4.6) about these fixed points gives three eigenvalues 



0, 3(2-7), 3 - y|'(A^x^ + A^x^) . 



(2.26) 



Thus the kinetic dominated solutions are unstable. 



F : The fluid dominated solution also exists for any A<^ and A^. There is no dark energy 
contribution. The eigenvalues for the system are 

3 3 3 , 
-2' -2' 2' ^'-''^ 
which show that the solution is unstable. 

SP : The dark energy dominated solution exists for A^ < y^A^ + 6. This means that as long 
as the slope of the potential of a scalar field X^f, is not much stiffcr than that of LW field A^ we 
can have the dark energy dominated solution. The eigenvalues for the system at this point are 

A^ A^ 

-3+^, -3 + ^, \lff-^l, (2.28) 
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which indicate that the solution is stable when the constraint A^^j < 87 is satisfied. Unlike a 
single scalar field models, one can obtain the accelerated expansion even if the potential is not 
too flat. To obtain the accelerating expansion we need that the slope of a potential of a scalar 
field should not be too much flatter than that of a LW field. Similar idea is adapted in so called 

"assisted quintcssece" [25, 26]. 

SPF : The dark energy fluid dominated solution exits for X^jjr > 87. This solution tracks the 
background barotropic fluid. Thus wc need a stiff scalar field potential compared to that of LW 
field. At this point, the eigenvalues for the system are 



-^(2-7), -^(2-7) 



1 + 



87(37-A|^ 



-4(2-7) 



1 I 87(37 -A^//) ' 
\ (2-7)A^,, . 

(2.29) 

Thus, we have the stable solution as long as the constraint 87 < A^^j < 247^/(97 — 2) is satisfied. 



3 Evolution without Self Interaction 

In the previous section, wc study the critical points and their existences and stabilities. The 
effective eos of the dark energy is given by wde = — 1 + X^j^/S at the point SP. Thus, we should 
have Ag^j < 2 in order to get the accelerating expansion at late time {i.e. o^de < —1/3) when 
the dark energy dominates the Universe. We can see the situation more clearly if we express the 
effective dark energy energy contrast Ode and its eos wde by using the dimensionless variables in 
(2.18) and (2.19) 



^^DE 



^2 



,,2 



x\ + 



(8.1) 
(8.2) 



Prom these equation, we can find the condition for late time accelerating universe. Energy density 
of dark energy component should be dominant and its eos should be smaller than —1/3 to get 
the late time accelerated expansion. In Table 3, we distinguish a scalar field dominated case 
(w^^ > — 1) with a LW field dominated cases (wj^^ < —1) for the specific choices of A^ and A,^. 
We mean the present values of any quantities by using (0) in the affix. For example, we obtain 



> x'^' from the above equation (8.2) to get a scalar field dominated universe at present. This 
condition can be satisfied if the slope of the potential of a scalar field lA^I is stiffer than that of 
a LW field |A^|. In this case a scalar field rolls a potential faster than a LW field and satisfy the 
above condition. Late time phantom dominated epoch wde < — 1 also can be reached when a 
LW field moves faster than a scalar field. This case will be happened when |A^| > lA(j,|. Only the 
magnitude of slope is concerned in those cases. Wc have the opposite signs of the slopes for scalar 
fields in LW theory as shown in the equation (2.8). When [A^l = |A^|, the scalar field (p and a LW 



„(0) 
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Table 3: The present value of dark energy eos wde for the various conditions for A<^ and X^. The 
slopes of potential determine the values of and at present giving the different values of 0;^^^. 
In this table and x^ represent the present values of them. 



'^DE 





< -1 


-1 


> -1 






< |Av.| 




> |A^| 






<4 




>< 



field t/j evolves same as each other. In this case, we have x'^ = x"^ and the eos of dark energy is 
always satisfied wde = — 1 even though each scalar field evolves. 

We show the evolutions of energy density contrast of dark energy, that of the matter, and the 
eos of dark energy for the specific values of and . 

Case I ) |A<^| > |A^| : 

In the left panel of Fig. 1 we show the evolutions of dimensionless variables x^, x^, and y 
when lA^I = 1.0 > |A^| = 0.8. In this case, a scalar field rolls down a potential faster than a 
LW field. Thus, we have bigger x^ {i.e. (j)') value than x^ {i.e. as shown in the left panel of 
Fig. 1. The y value is also increased as fields evolve. We can understand this evolutions from 
equations (2.22) - (4.6). All of x^, x^, and y values are increased as they evolve in the universe. 
In the X(j)-x^ plane, we can see that x^ starts from and reaches to 0.41 while x^ increases from 
to 0.33. In both x^-y and x^-y planes, we find that y is increased from to 0.97. Because x^ 
is bigger than x^/,, y is smaller than 1 as can be seen in the equation (3.1). Also the evolutions 
of the energy density contrast of dark energy iloE, that of dark matter ilm, and the eos of the 
dark energy wde are shown in the right panel of Fig. 1 for the same values of A^ and A,^ . As we 
show in Table 3, in this case the present value of the dark energy eos is bigger than —1. The eos 
of dark energy evolves starts from and reach to —0.88 at present. Also the energy density of 
dark energy takes over that of matter as they evolves. When we choose the positive value of A^, 
the above argument is same except the sign of x^. In that case, x^ decreases as it evolves with 
the same absolute values as the negative value of A^ case. Thus, the cosmological values like J^de 
and wde will not be changed by changing the sign of A^ . 



Case II) |A^| < |A^| : 

We also show the evolutions of x^, x^, and y for lA^I = 0.8 < |A^| = 1.0 in the left panel of 
Fig. 2. In this case, a scalar field climbs a potential slower than a LW field. Thus, x^ value is 
smaller than x^ value. Again y- value is increased as fields evolve but it can reach to bigger than 1 
at present. This happens because in this case x^ < x^ and in order to get the Ode < 1 we should 
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Figure 1: a) The evolution of x^, x^, and y for A,^ = 1.0 and = —0.8. b) The evolution of 
^^DE) ^m, and LOBE for the same values of A,^ and A^ . 

have y > 1. In the x^-x^ plane, we see that both x^p and x^ increase from to 0.33 and from to 
0.41, respectively. Also from both x^p-y and x^-y planes, we see that y starts from and reaches 
to 1.03. The evolutions of f^DE) ^m, and lude are indicated in the right panel of that figure for 
the same values of A^ and A^. In this case, present value of dark energy eos is smaller than — 1 
as we can find in Table 3. The dark energy eos evolves from and reaches to —1.12 at present. 
Whenever, the condition lA^I < |A^| is satisfied, we have phantom attractors at late time. 



In both Fig 1 and 2, we show the evolution of a scalar field and a LW field in the specific 
values of A^ and A^. However, the present energy density of dark energy component is too big 
in both case {^y)E — ^) compared to the observed value (^^^e — 0.73). We can see this situation 
from equations (3.1) and (3.2). In order to get the dark energy eos close to —1 we obtain ~ 
from the equation (3.2). If we use this condition to the equation (3.1), we get Ode = y^- We 
also have V ~ V^_^ from the equation (2.4) because we have A^c/) = — A^V in this case if (p and 
ip started from same initial values. Thus, JIde — 1 in this case and can not be matched to the 
observed value Ode — 0.73. 

We need to tune the parameters to obtain the suitable current Ode value. In order to get 
^de ^ should have A,^ > A^ as shown in the above. As an example we can have Ode = 0.73 
when we choose A,^ = \/4.5 and A^ = VOA as shown in Fig 3. This choice of values satisfies the 
existence of dark energy-fluid dominated solution as shown in Table 1. Nevertheless, this case 
follows the matter tracking attractor and the eos of dark energy converges to 0. Thus, this case 
will not be able to derive the late time accelerating universe even though it dominates the late 
time energy component. 



-8 -6 -4 -2 

In a 



Figure 2: a) The evolution of x^, x^, and y for = 0.8 and = —1.0. b) The evolution of 
r^DE, ^m, and wde with the same values of A^ and A^ . 



y 




In a 



Figure 3: a) The evolution of x^, x^, and y for A,^ = \/45 and A^ = — \/0~4. b) The evolution of 
f^DE) ^mi and cjde for the same values of A,^ and A^ . 
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Ode 








(ODE 



-8 -6 



■A 

Ina 



Figure 4: a) The evolution of I^dEi and wde for A,^ = 4.1, = —0.4, and M = 10 Hq. b) 
The evolution of same quantities for Xa, = 2.1, A^, = —0.1, and M = W~^Hq. 



4 Evolution with Self Interaction 



In this section, we consider the model with including the self interaction term of LW field, 
M'^ijP' /2 as in the equation (2.3). As in the previous section, we introduce new dimensionless 
variables similar to equations (2.18 - 2.19) including the self interaction term 



Kip' HiM%l) 

. = (4.3) 

We can find the evolution equations of the above variables (4.1 - 4.3) with using the interaction 
term as given before (2.4) 



H' 



= 


-(if- 


y' = 


[H' 


x'tp = 





M 



(4.4) 
(4.5) 
(4.6) 
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Compared to the without self interaction case, we are not able to have the critical points in this 
system. However, we can still investigate the cosmological evolution of this system. 

In the left panel of Fig 4, we show the evolution of OdEj and cjde when we choose that 
= 4.1, = -0.4, and M = IQ^^-'Hq. In this case, we have = 0.73 and = -1.83. 
Thus, this case can be viable for the current observations. In the right panel of this figure, we 
show the cosmological evolutions of same quantities for A,^ = 2.1, A^ = —0.1, and M = 10~^i/o- 
With this specific choice of parameters we have Oj^^ = 0.73 and uj^y^ = —0.01. This is not be able 
to produce the late time acceleration universe even though the dark energy component becomes 
a dominate component in the universe at late time. Also as we can see in the figure, the energy 
density of the dark energy dominates the universe too early. 

5 Conclusion 

Wc show a dark energy candidate model from LW theory. With two scalar fields we can have the 
late time equation of state of the dark energy across —1 which might be suitable for the current 
observations. Even though the structure of the theory is similar to the so called "quintom" model, 
wc can not avoid the interaction between two scalar fields when wc start from LW theory. Also 
the signs of slope of potential are opposite to each other. Interestingly, there exist the late time 
tracking attractor solutions for the specific choice of the slopes of a potential. However, this simple 
exponential potential models can not give the viable models which derive the current late time 
accelerating universe. We also investigate the model with including the self interaction term. In 
this case we can produce the suitable cosmological evolution of observable quantities. However, 
we are not able to find the stable solution in this case and we also suffer from the fine tuning for 
the mass of LW field. 
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